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ABSTRACT 


An experimental investigation was made of the modulation of an 
acoustic standing-wave field due to trapped, cavitation-induced bubbles. 
Maximum acoustic pressures generated are on the order of 0.1 bar. It 
is verified that bubbles driven below resonance are trapped just above 
pressure antinodes of the standing-wave field; bubbles driven above reso- 
nance, just above nodes. Maximum modulation percentages of the standing- 
wave field by bubble scattering are found to be in the neighborhood of up 
to 12%. Modulation noise spectrum levels are recorded and are found to 
decrease from 2.5 Hz to 100 Hz with slopes ranging from -3.2 to -5.4 dB/ 
octave. Less than 1% of the modulation noise energy lies in the region 


above 100 Hz. 
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I. IDNTRODUCTION 


This paper concerns the effects of cavitation-induced bubbles in ampli- 
tude modulating a medium intensity (from about 0.02 bar (ms) to about 0.15 
bar (xms)) standing wave field. The bubbles, after generation, are trapped 
at the nodes or anti-nodes of an acoustic standing wave field, and at suffi- 
ciently high driving levels dance around. Because all of the bubbles gene- 
rated are grouped fairly closely in size around that value which would be 
resonant at the exciting field frequency, they act as good sound scatterers. 
The constant generation, growth in size, dancing about and eventual escape 
of the trapped bubbles causes the scattered acoustic energy to modulate 
randomly the exciting carrier signal. The purpose of this research was to 
measure the frequency spectra for this random modulation signal. 

Preliminary to development of these spectra, however, it was necessary 
to investigate the nature of the standing wave field, which led directly 
into another interesting and as yet little explored area, namely, the nature 
of individual bubble motion or dynamics. It was felt that a description 
(largely qualitative due to equipment and time limitations) of individual 
and aggregate bubble behavior would be useful, both as an area for study in 
itself and as providing a possible clue to the nature of the observed sound 
eeectra. 

Thus, the investigation resolved itself into three broad areas, nanely, 
the nature of the standing wave field, the development of the scattered 
noise modulation spectra, and the nature and characteristics of bubble motion 


within the field. 
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Although quite a bit of research had been conducted on cavitation 
bubbles and cavitation noise, the above areas have been relatively 
unexplored. In some early observations, Boyle (Ref. 1) and Hopwood 
(Ref. 2) demonstrated that bubbles tend to collect at the pressure nodes 
of a vertical standing wave pattern. Samewhat later, Blake (Ref. 3) 
showed that bubbles would also collect at pressure anti-nodal positions, 
provided that the resonance frequency of the bubble is greater than that 
of the applied sound field. Eller (Ref. 4) provided a theoretical 
explanation for the trapping process and calculated the minimum pressure 
amplitude required to trap a bubble. Strasberg and Benjamin (Ref. 5) 
observed that if the applied field were strong enough, the trapped 


bubbles would dance about erratically. 
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II. THEORY 


This section contains a theoretical discussion of the standing wave 
field and of certain aspects of bubble dynamics. A rigorous treatment 
of the exceedingly canplex modal structure within the tank and of the 
minute erratic motions of the visible bubbles under the influence of the 
acoustic field is beyond the scope of this paper. Therefore, the following 
assumptions are made: 

1) The presence of pressure release boundaries at all surfaces. This 
is true at the water surface, but only approximately so at the sides and 
bottom of the tank, which are of thin glass. 

2) Perfect perpendicularity of the tank sides, and that the water 
surface is parallel to the tank bottom, neither of which is exactly true. 
3) The non-existence of any discrete lines in the spectrum of the 

noise modulation signal. In fact, none were detected; however, weak 


frequency camponents may have been present. 


A. THE STANDING WAVE FIELD 

The following discussion follows that in Morse (Ref. 6). We consider 
a rectangular tank measuring 60 cm by 35 cm in length and width respec- 
tively, and filled to a depth of 17 cm with water. A rectangular coordinate 


system is established as in Figure 1 below. 
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The wave equation in rectangular coordinates is 
2 2 Z 
d*p , 24p , ap _ _1 3*p (1) 
ax ay? = az? ce# 82 
where x, y, and z are spatial coordinates, and p is acoustic pressure. 


Solutions to Equation (1) are of the form 


= bad : Ne So ons 
Pp = p.sin(o,— + Se i $,) Sin (we + $,)exp(-2nift) (2) 
where f = 5 lw, + Hy, + w, ] and dt ty and o, are determined by the 


boundary conditions. If we assume that all boumdaries are pressure 


release, the o's are all zero and we have 


= j io } ve , cS = , 
P = p,Sintu = dsintus- )sin(u,clexp (-20ift) (3) 
ec _ me Bee, 
where w= xT m,, = lj2;o--..-. oltnilarly, we have Wy = yy" 
Ba = 
Ww, = ze" where M9 Mm yoo accwom Loe Mes M+ m, are related by 
ne 2 1 3 1 ee, 
FF + +S) (4) 
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Each of these solutions describes a free vibration mode of the enclosure. 
Thus the wavelength measured along any axis must at all times be greater 


than on as can be demonstrated by substituting values for mo Mm, and mM, 
A 
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and calculating the equivalent yo A. and an The lowest frequency mode 


y 
allowed in the tank corresponds to the (111) mode, where the triplet of 
numbers represent gH, tbe respectively. Applying (4), we have the 
result that this mode corresponds to a frequency of 4.18 KHz. Below 
this frequency there can exist no standing wave structure in the tank. 
At discrete frequencies above this cutoff, resonant modes can exist, 
these frequencies being determined by the eeoulere a oh AL Thus the 
first several frequencies above the (111) mode, in ascending order, are 
those corresponding to the triplets (211), (311), (121), (221), (321), 
Gio), (411), etc. 

At higher frequencies, the number of modes that can exist within a 
given frequency band becomes greater and greater. If we treat the Wt 
ws Ww as camponents of a three dimensional vector in frequency space it 
can be shown that the number of modes dN contained within a frequency band 


af is given by 
a Anf£2v - 


c3 


oN (5) 


where V = a Yo 2 is the volume of water in the tank, 
so that for our particular tank of water being ensonified with a 30 KHz 
Sinusoidal signal, we should have dN = 12 for a df of 100 Hz. Thus, 
within a 0.3% frequency spread around 30 KHz we should have about 12 
possible modes of vibration. 

The driving transducer forces vibrations in the tank, with the result 
that all modes of vibration can be excited, but those modes lying closest 
to the exciting frequency are more strongly excited than those further 


from it. The relative amount of energy present in any mode depends upon 


the Q of the system, where Q is defined by 








where £ is the exciting frequency and 

f, and f. are the lower and upper frequencies, respectively, at which 
the acoustic pressure amplitude is down by 3 dB. A system such as our 
tank, with a relatively high Q (Q = 200-500), will have most of the 
vibrational energy concentrated in only a few of the allowed modes 
grouped around the exciting frequency. The lower the Q of the system, 
on the other hand, the more evenly will the energy be distributed over 


a larger number of allowed modes. 


B. BUBBLE DYNAMICS 

A bubble located in a standing wave acoustic field will experience 
a net force from that field. This force is due to spatial pressure vari- 
ations within the field. If bubble pulsation and acoustic pressure are 
out of phase, that is, if the bubble grows smaller as acoustic pressure 
increases, then the bubble will tend to be driven toward a pressure 
antinode. If the reverse condition holds, with bubble size and acoustic 
pressure being in phase, the bubble will tend to be driven toward a 
pressure node. These phenomena may be utilized to trap bubbles at 
relatively stationary positions in an acoustic standing wave field. It 
can be shown (Ref. 7) that the minimm acoustic pressure amplitude at an 
antinode necessary to trap a bubble in a standing wave field in the z 


direction is given by 


2 (1-22)2 + d2 
SI guid ll®) So Cl Caoreecerenmereica (7) 
nf ees) 
= pw*Ro* and 1S equal to 1.0 at resonance, 


2 
where & 3p, 


R, is the equilibrium radius of the bubble, 
d is a damping term dependent on the contributions due to heat 
Gonduction, sound radiation, and viscosity, 


9 








n is a number between 1.0 and 1.4, depending on whether the bubble 
oscillation process is essentially isothermal or adiabatic, 
and Poin refers to the pressure amplitude (Zero-to-Peak) at a pressure 
antinode of the wave field. 

In the above result, we have assumed that the wavelength in the z direction 
has its smallest possible value, re Then (7) gives us the minimum pres- 
sure at an antinode at which a bubble will be trapped in the acoustic field 
against the force of buoyancy. Figures 2 and 3 give minimum trapping pres- 
sure vs bubble size for some typical frequencies. The significance of these 
figures is to fix a rough minimm value of the pressure amplitudes at the 
anti-modes in the system during the experiment. At a given frequency, the 
acoustic pressure necessary to trap a bubble has minimum values just above 
and below resonance, and increases as we move away from resonance in either 
direction. Bubbles driven below resonance will be trapped, relative to the 
Z axis, Just above pressure maxima, while bubbles driven above resonance 
will be trapped just above pressure nodes. 

Regarding the trapping process, it is interesting to note that if a pure 
modal state or cambination of modal states could exist within the tank, there 
could be no trapping at frequencies higher than resonance, since modal posi- 
tions within the tank would be vertical and horizontal intersecting planes, 
and bubbles, once driven to a vertical plane, would experience no further 
vertical acoustic force, and would escape. However, due probably to the 
absorbtion of acoustic energy at the walls and various perturbation effects, 
such as the presence of the transducers (the driving transducer measured about 
3 cm in diameter and 4 am in length), and non-parallelness of the boundaries, 
the theoretically predicted vertical nodal planes were not detected in the 
sound field, about which more will be said later. ‘Thus, bubbles could be 


trapped above resonance. 
10 
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IIIT. EXPERIMENTAL PROCEDURE 


The arrangement of equipment was campletely straightforward and is 
illustrated in Figure 4. The sine wave from the signal generator was 
amplified and used to drive a PZT-4 transducer. A variable inductor was 
connected in series to balance the capacitance of the transducer and 
provide a more effective impedance match with that of the amplifier. 

The acoustic signal was received on a small 1/8" diameter cylindrical 
probe hydrophone, given 40 dB of amplification, passed through a 60 Hz 
filter, then either viewed directly or passed through a rectifier and 
detector and fed to a 1/3-octave frequency analyzer. Late in the experi- 
ment it was discovered that a significant amount of the carrier was being 
leaked through the detector, so an additional low-pass filter, a Kronh-Heit 
active type, was installed after the detector. 

The transducers were installed in a standard 15 gallon fishtank, 
measuring 60 x 35 x 30 cm in length, width, and height respectively, and 
filled to a depth of 17 cm with water. The driver was clamped permanently 
in one corner of the tank and not touched through the experiment, while 
the receiving transducer was mounted on a Sliding carriage that permitted 


it to cover most of the tank, at least horizontally, with facility. 


A. MEASUREMENT OF NOISE 

The receiving hydrophone was calibrated at the start of the experiment 
by direct comparison with a calibrated laboratory standard. All pressure 
amplitudes calculated during the course of the experiment were figured 
directly from the results of this calibration. At the close of the experi- 


ment, the receiving hydrophone was recalibrated. The results of the two 
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calibrations were within 20%. In general, the pressure measurements con- 
tain the largest percentage errors of any of the quantities calculated 
during the experiment. Their accuracy (pressure amplitude measurements) 
1s estimated at + 203%. 

A typical unfiltered signal, taken directly from the hydrophone (after 
40 dB of amplification) is illustrated in Figure 5. Our main interest 
lay not in the carrier but rather in the irregular amplitude modulation 
of the carrier. 

After the removal of any 60 Hz component present, the signal was 
passed through a detector, basically a low-pass filter and rectifier, with 
a -3 dB point of about 1200 Hz. What remains is the amplitude modulation 
signal as pictured in Figure 6. 

Very little could be conjectured about the frequency spectrum of the 
amplitude modulation signal before it was analyzed, so it was fed intoa 
GR 119 wave analyzer, an instrument with a practical low frequency cutoff 
of about 20 Hz. It was soon apparent that by far the greatest percentage 
of the energy in the noise signal lay below that frequency. Therefore, 
the GR 119 was replaced with a GR 1564 octave band analyzer, which has a 
low-frequency cutoff of 2.5 Hz, and the noise measurements were made on 
that. As opposed to the GR 119, which had a constant width gate (3,10 or 
50 Hz), the GR 1564 was capable of operating in a 1/10 or 1/3 octave mode; 
the latter mode was used in the experiment. However, the GR 119 was used 
during several runs (using the 3 Hz gate) to search for discrete components 
in the noise signal. 

The 1/3 octave band level readings obtained were corrected to spectrum 
level using the conversion: 


SPHCTRUM LEVEL (dB) = BAND LEVEL (dB) ~ 10 log Af 
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where Af is the frequency spread of the band in question. Spectrum level 
is then the average amount of energy per hertz bandwidth that one measures 
in any band. The GR handbook has a table of Af values for the 1/3 octave 


band used in the experiment. 


B. OBSERVATION OF BUBBLE MOTION 

The original experiment was not set up for the purpose of investigating 
bubble dynamics within the tank, and the apparatus and general setup did 
not lend themselves to the taking of precise, quantitative data in this area. 
However, the bubble dynamics proved so interesting in itself that it was 
decided to include a section describing, qualitatively, bubble effects, and 
some rough calculations of bubble size as related to trapping pressures. 

The bubble generation and trapping process was observed through the tank, 
Side illuminated against a black backcloth. This system provided excellent 
contrast and good viewing. 

A bubbie's size was determined from its velocity of ascent through the 
water by using Stokes' Law. The assumption was made that the bubbles' 
shapes remained spherical as they rose, and no attempt was made to apply 
empirical corrections to Stokes' Law calculations, since these corrections 
were small compared to the error in determining the bubble's position and 
its time of rise. The procedure used was to let the bubble field stabilize 
for a few minutes, then to select a bubble, and by holding a ruler externally 
against the tank, determine its vertical position. Parallax errors were 
reduced to as small a value as possible by eyeball. The sound field was 
then rapidly shut off; a stopwatch being started at the same time. The 
bubble's time of ascent to the surface was then determined, the velocity 
thereby calculated, and Stokes' Law, in the following form, employed to find 


bubble radius: 
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R= (52) (8) 


where v, the kinematic viscosity, = 0.0093 poise, appropriate for 23°C 


tI 


g = 980 an/sec 
u = bubble velocity 
R = bubble radius in cm 


A plot of Rvs u is given in Figure 7 below. 
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12 STOKES' LAW PREDICTION OF 
BUBBLE SIZE VS VELOCITY 
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C. CALCULATION OF Q 

Typical values for Q were calculated by tuning the frequency generator 
to get a relative peak in the carrier signal, then fine tuning it to locate 
the -3 dB points. This procedure is valid only as long as a single mode is 


present over the range of frequencies in question; that is, from the lower 
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3 dB point to the higher 3 dB point. Typical values for this frequency 
spread ranged from 50 Hz to 250 Hz. Refer to Table 1. Our previous cal- 
culations, however, indicate that at 30 KHz, we should expect about 12 
vibrational modes per 100 Hz bandspread. In actuality, although some 
modes were extremely close together, too close to permit a calculation Q, 
there were numerous modes, which, to the eye, appeared to be the only 
Significant mode present over a bhandspread of several hundred hertz, and 
thus pemnitted calculating values for Q. It is not understood why so 
many fo the theoretically possible model configurations are absent, or at 


least, unmeasureble. 


D. MEASUREMENT OF MODULATION PERCENTAGE 

Modulation percentage, as defined here, is the ratio of the rms fil- 
tered modulation signal to the rms modulated carrier. The mms filtered 
modulation voltage was measured on the GR-1564 octave band analyzer, using 
the all-pass band. During measurements of modulation percentage, the 
hydrophone was located at a position in the sound field where the carrier 


was a relative maximum. 


E. MEASUREMENT OF PRESSURE AMPLITUDE 

The receiving hydrophone was calibrated by comparison with a laboratory 
Standard. The receiving Signal, regardless of further processing, was 
given 40 dB of amplification. The actual receiving hydrophone sensitivity 
was 0.01 volt (rms) = 0.016 bar (rms). However, all scales on the graphs 
included herein represent actual meter readings, that is, with the addition 
of 40 GB amplification. Therefore, we have that 1.0 volt (xms) = 0.016 


bar (mms). 
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IV. EXPERIMENTAL RESULTS 


A. ACOUSTIC WAVE FIELD 

The sound field, in spatial configuration, when probed along a line 
parallel to one of the boundaries at frequencies above the low-frequency 
cutoff, demonstrated the expected sinusoidal pattern. At relatively low 
frequencies the mode shape came close to agreeing with the theoretical 
description of a pure mode. At higher frequencies however, the mode shape 
appeared to be the sum of several modes of oscillation. Typical field 
patterns are shown in Figure 8 and 9. These patterns were made at low 
driving levels so that the effect of modulation noise was minimal. The 
expected sinusoidal shape is evident, although the differing heights of 
the various peaks indicate the presence of more than one mode in the tank. 
Note also that at no node does the acoustic pressure amplitude go to zero, 
indicating that the walls of the tank are not true non-absorptive pressure 
release boundaries. 

An attempt was made to construct the (111) mode within the tank, 
without success. The lowest obtainable modal configuration was the (311) 
mode, at a measured 5.3 KHz. 

It is difficult to give representative values of Q. As the frequency 
generator was swept over a range of frequencies from, say, 10KHz through 
40 KHz, with the receiving hydrophone placed at an arbitrary spot in the 
tank, there were observed hundreds of stronger and weaker resonances. Many 
of these resonances were so close to each other that the separation between 
them was less than the frequency spread necessary to measure the -3 dB 


points. Several resonances at various frequencies are tabulated in 
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Table 1 below. The three values at each frequency correspond to 
different positions in the tank. For a given mode, then, the values for 


Q do not differ by more than about 20%. 


TABLE 1: Representative Values of Q 


Frequency (KHz) Af Q. 
25.711 54 474 
69 372 

60 429 

27.050 120 226 
105 258 

136 198 

30.207 163 197 
170 178 

191 158 

sul, yal 130 240 
132 254 

144 218 

32.543 102 319 
115 282 

al 269 


In the theoretical calculations it was determined that there should 
exist about 12 possible modes within a 100 Hz band centered at 30 KHz. In 
fact, only two or three are observable. Changing the positions of the 
sending or receiving hydrophones <did not increase the number of observed 


modes at a given frequency. 


B. MODULATION SPECTRUM 

Figure 6 shows a typical output signal, taken directly from the 
receiving hydrophone, and given 40 dB amplification. Here V. = 70 volts 
(xms) which corresponds to the maximum driving levels used during the 


experiment. Cavitation was severe and the bubble field was stabilized. 
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The low frequency modulation is apparent. Figure 10 shows a representa- 
tive unfiltered carrier frequency spectrum, with 40 dB of amplification. 
Here we are also driving at a high level, and the modulation percentage 
is about 10%. The lower hum in Figure 10 is a measure of the 60 Hz 
noise interference in the output signal from the amplifier, while the 
upper hump is the carrier signal and its noise modulation products on 
either side. 

Measurements of the modulation percentage aReote that at the maxi- 
mum driving levels used, the modulation percentage was on the order of 
8-12%. The filtered modulation signal corresponding to the unfiltered 
Signal in Figure 6 is shown in Figure 7. The vertical sensitivity of 
the CRO has been increased to better delineate the sicnal. Reference to 
the horizontal sweep time rate indicated that a large amount of the 
Signal's energy lies at very low frequencies (less than 10 Hz). 

Typical noise modulation spectra are shown in Figures 11, 12, and 13. 
As can be seen, the spectrum level in every case diminishes in magnitude 
from the lower to the higher frequencies. In no case did the amount of 
energy in the modulation spectrum above 80 Hz amount to as much as 13% of 
the total energy in the rectified signal. One apparent discrepancy evi- 
dent in Figures ll, 12, and 13 is the fact that the total energy of any 
Signal measured was 2-3 dB greater than the sum of the individual bands 
added together. The only explanation offered is that this could be an 
idiosyncrasy of the GR 1564 or might be a consistent error in data 
recording, namely, a constant lower-than-actual reading of the dB meter. 
This could amount to a considerable error in the summation of 17 bands, 


which was the number used. During measurement of these modulation spectra, 
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the maximum carrier mms pressure amplitudes were typically in the order 
of up to 0.04 bars. The percentage modulation was typically from 4% to 


12% at positions of maximum carrier amplitude. 


C. BUBBLE DYNAMICS 

As mentioned earlier, bubble trapping occurs at positions just above 
either vertical nodes or antinodes, depending on whether the bubble reso- 
nance frequency is lower or higher than the acoustic exciting field fre- 
quency. ‘The lowest frequency at which trapping occured was 5.8 KHz. At 
this frequency, the vertical pressure amplitude profile consisted of nodes 
at the surface and bottom of the tank and an antinode in the middle. Thus, 
only bubbles below resonance were trapped (just above an antinode). At 
sufficiently high driving levels, the bubbles began to grow (by rectified 
diffusion) and finally escaped to the surface. The time fran trapping to 
final escape, at this frequency, was on the order of 15-30 seconds at the 
highest driving levels used. 

At the driving frequencies used for most of the study, that is, above 
about 25 KHz, the vertical pressure pattern was no longer a single half- 
wave, but several half-waves, thus permitting nodal as well as anti~nodal 
trapping. In fact, the higher the driving frequency, the shorter the time 
time period a given bubble would stay small enough to be trapped in an 
anti-nodal position. At frequencies above 30 KHz, the time pericd between 
observed bubble generation (due to cavitation) and the jump from an anti- 
nodal to a nodal trapped position was not more than several tenths of a 
second and could not be quantitatively measured with the available equipment. 
Bubbles trapped above their resonance frequency (that is, at nodal positions), 
were generally stable, and depending on the driving voltage amplitude, would 


maintain their positions for several minutes at a time. 
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The number of trapped bubbles in the tank at a given time was depen- 
dent upon the driving frequency, driving voltage amplitude, and the length 
of time the field had been activated. At an estimate, the number varied 
between 15 and 50 bubbles. Generally, there seemed to be no preferred 
nodal or antinodal positions, and the trapping process occured fairly uni- 
formly throughout the tank. Quantitatively, the all-pass noise modulation 
Signal was measured for several different frequencies and driving levels, 
as a function of receiving probe position. Under no combination of driving 
frequency and amplitude did the variation within the tank exceed 2 dB. 

Given a trapped bubble at a fixed frequency, the increase in its motion 
at higher driving levels was quite marked. At lower levels, the trapped 
bubbles were almost stationary, merely undergoing small excursions about an 
equilibrium position. As the driving voltage was increased, the motion 
became more violent, following an ellipsoidal path. At the highest driving 
levels obtainable, the bubble motion amounted to excursions of as much as 
a centimeter (estimated by eye). At these levels, the acoustic pressure 


amplitude amounted to on the order of 0.15 bar (rms). 


D. MEASURED BUBBLE SIZES 

Bubble sizes calculated through application of Stokes’ Law revealed that 
the smallest measured, trapped bubble had a radius of 85 microns, in an 
exciting field of 30.744 KHz. Resonance size at this frequency was 104 
microns. Thus, at a given frequency, there was a large variation in bubble 


emze on each side of resonance. 
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V. CONCLUSIONS 


1. Bubbles trapped in a stationary acoustic field cause scattering of 
acoustic energy and can amplitude modulate the ensonifying acoustic field. 
2. Bubbles can be trapped below and above resonance, near anti-nodes 
and nodes of the standing wave field, respectively. 

3. Carrier signal pressure amplitudes of about 0.07 bar (rms) produced 
severe bubble motion. 

4. At pressure amplitudes ranging up to about .04 bar (mms), the 
amplitude modulation percentage was about 8-12% when bubbles were present. 
5. For the frequencies and pressures observed, the modulation spectra 
levels due to the scattered noise decreased with increasing frequency, 
with slopes that ranged from 3.2 dB to 5.4 dB per octave, over a range 


of frequencies from 2.5 Hz to 100 Hz. 
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